Actomyosin contractility rotates the cell nucleus 
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The nucleus of the eukaryotic cell functions amidst active cytoskeletal filaments, but its response 
to the stresses carried by these filaments is largely unexplored. We report here the results of studies 
of the translational and rotational dynamics of the nuclei of single fibroblast cells, with the effects of 
cell migration suppressed by plating onto fibronectin-coated micro-fabricated patterns. Patterns of 
the same area but different shapes and/or aspect ratio were used to study the effect of cell geometry 
on the dynamics. On circles, squares and equilateral triangles, the nucleus undergoes persistent 
rotational motion, while on high-aspect-ratio rectangles of the same area it moves only back and 
forth. The circle and the triangle showed respectively the largest and the smallest angular speed. 
We show that our observations can be understood through a hydrodynamic approach in which the 
nucleus is treated as a highly viscous inclusion residing in a less viscous fluid of orientable filaments 
endowed with active stresses. Lowering actin contractility selectively by introducing blebbistatin 
at low concentrations drastically reduced the speed and persistence time of the angular motion 
of the nucleus. Time-lapse imaging of actin revealed a correlated hydrodynamic flow around the 
nucleus, with profile and magnitude consistent with the results of our theoretical approach. Coherent 
intracellular flows and consequent nuclear rotation thus appear to be a generic property that cells 
must balance by specific mechanisms in order to maintain nuclear homeostasis. 

Keywords: Nuclear dynamics — Cell shape — Actomyosin — Rheology — Hydrodynamics of active gels — 
Active stresses 



INTRODUCTION 

The nucleus is the largest and stiffest organelle in a eukaryotic cell [1]. It is actively coupled to the dynamic 
cytoskeleton [2-5] by means of a variety of scaffold proteins: contractile [6] acto- myosin complexes, microtubule 
filaments constantly undergoing dynamic reorganization, and load bearing intermediate filaments [2-4, 7-10]. The 
nucleus has been found to translate and rotate during cell migration [11—17]. It is reasonable to suppose that such 
motions are a result of active processes in the cytoplasm, involving the cytoskeleton and molecular motors [2, 3, 7- 
10, 12, 14, 17]. The positioning of the nucleus in the cellular environment is critical to many physiological functions 
such as migration, mitosis, polarization, wound healing, fertilization and cell differentiation [16, 18]. Alterations to 
nuclear position have been implicated in a number of diseases [16, 18]. Taken together these studies suggest that 
the mechanical homcostatic balance of nuclear positioning and dynamics is intimately coupled with cellular geometry. 
While a number of molecular players have been implicated in this context [2-5, 9, 10], the role of actomyosin 
contractililty on nuclear dynamics has not been explored. 

In this paper, we show that cell geometry and active stresses are critical components in determining nuclear 
position and movements. Fibroblast cells (NIH3T3) plated on micro- patterned fibronectin surfaces of varying shapes 
and aspect ratio were used to assess the effect of geometrical constraint on the translational and rotational movement 
of the nucleus. Time-lapse imaging revealed a correlation between actin flow patterns and nuclear movement. We 
show that a hydrodynamic model of oriented filaments endowed with active contractile stresses [19-23], with the 
nucleus entering only as a passive inclusion, gives rise to the observed organized actin flow and nuclear rotation. 
While preparing the present work for submission, we became aware of two works [24, 25] with theoretical formulation 
and predicted behaviours similar to ours. The contexts in which these works are set is different from ours, i.e., the 
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dynamics of the cell nucleus is not the subject of these papers. In addition, the boundary conditions are different 
in detail. Reference [24] was in a Taylor-Couette geometry, i.e., there is no medium inside the inner circle, and 
Reference [25] was in a circular geometry without a central inclusion. Our observations suggest that nuclear rotation 
and circulating flows are an inherent property of the active cell interior under geometric confinement. That nuclear 
rotation is not a normally observed feature of cell dynamics suggests that the cell must possess other mechanisms to 
suppress it. We discuss these towards the end of the paper. 

MATERIALS AND METHODS 

Cell Culture: NIH3T3 fibroblasts (ATCC) were cultured in low glucose DMEM (Invitrogen) supplemented with 
10% Fetal Bovine Serum (FBS) (Gibco, Invitrogen) and 1% Penicillin- Streptomycin (Invitrogen). Cells were main- 
tained at 37° C in incubator with 5% C02 in humidified condition. Cells were trypsined and seeded on fibronectin 
coated patterned surfaces for 3 hours before staining or imaging. For confocal imaging, low well ibidi non- treated 
hydrophobic dishes were used. 65,000 cells were seeded on each time on patterned surfaces (with 10,000 patterns) for 
30 minutes, after which the non-settled cells were removed and media was re-added in the dishes. Blebbistatin (Invit- 
rogen) were diluted from stock using filtered media. Blebbistatin was used at concentration of 1.25/iM. This minimizes 
the effect of any other solvent like DMSO. Microtubule was immunostained using a-tubulin antibody (1:200, Abeam) 
and Alexa fluor 546 secondary (1:500, Invitrogen) in cell plated on triangular pattern. The nucleus was labeled using 
Hoechst (1:1000). 

Preparation of PDMS Stamps and micro-contact printing: PDMS stamps were prepared from PDMS 
Elastomer (SYLGARD 184, DOW Corning) and the ratio of curer to precursor used was 1 : 10. The curer and 
precursor were mixed homogeneously before pouring onto the micropatterned silicon wafer. The mixture was degassed 
in the desiccator for at least 30 minutes to remove any trapped air bubbles and was then cured at 80°C for 2 hours, 
after which the stamps were peeled off from the silicon wafer. Micropatterned PDMS stamps were oxidized and 
sterilized under high power in Plasma Cleaner (Model PDC-002, Harrick Scientific Corp) for 4 minutes. 30/d of 
100^ig/ml fibronectin solution (prepared by mixing 27/d of lxPBS to 1.5/ig of lmg/ml fibronectin and 1.5/d of Alexa 
647 conjugated fibronectin) was allowed to adsorb onto the surface of each PDMS stamp under sterile condition for 
20 minutes before drying by tissue. The PDMS stamp was then deposited onto the surface of a low well non-treated 
hydrophobic dishes (Ibidi) (for high-resolution imaging) to allow transferring of the micro-features. Subsequently, the 
stamped dish was inspected under fluorescent microscope to verify the smooth transfer of fibronectin micro-patterns. 
Surface of sample was then treated with 1ml of 2mg/ml Pluronic F-127 for 2 hours to passivate non- fibronectin 
coated regions. 

Cell Transfection: Transfection of various plasmids in wt NIH3T3 cells was carried out using JetPRIME polyplus 
transfection kit. 1/xg of plasmid was mixed properly in 100/il of JetPRIME buffer by vortexing and spinning, 3.5/il of 
JetPRIME reagent was then added and the mixture was again vortexed and spun. The mixture was incubated for 30 
minutes and then added to 50-60% confluent culture in 35mm dish. Cells were kept in fresh media for 2 hours prior 
to addition of transfection mixture. Cells were incubated for 20 hrs before plating them on the patterned substrates. 

Imaging: Phase contrast imaging of cells on different geometrical patterns was done on Nikon Biostation IMq 
using 40x objective at 37°C in a humidified chamber with 5% CC>2- Confocal time lapse imaging of cells transfected 
with various plasmids (Lifeact EGFP, r RFP and dsRed ER) was carried out on Nikon AIR using 60x, 1.4 NA oil 
objective at 37°C in a humidified incubator with 5% CC>2- 

Image Analysis and quantifications: Acquired images were processed and analysed using ImageJ software 
(http://rsbwcb.nih.gov/ij/index.html). To determine the translational coordinates and rotational angle of nucleus, 
diagonally opposite nucleoli were manually tracked from the phase contrast image of the cell using the ImageJ plugin- 
MtrackJ (http://www.imagescience.org/meijering/software/mtrackj/). The translational and rotational autocorrela- 
tion were calculated from the residual of the linear fit to corresponding curves- the detrended curves, thereby taking 
into account only the time scales relevant in our measurements. Particle image velocimetry (PIV) analysis was carried 
out using Matlab PIV toolbox-Matpiv between consecutive image frames separated by lmin. Images acquired were 
512 X 512 pixels. The size of the interrogation window was chosen to be 32 X 32 with an overlap of 50% between the 
consecutive time frames. The single pass method was used for calculating the velocities. Quantifications were done 
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using custom written program in either Lab VIEW 6.1 or MATLAB R2010a. All the graphs and curve fittings were 
carried out using OriginPro 8.1 (OriginLab Corporation, Northampton, USA). 



RESULTS AND DISCUSSION 
Geometric constraints on the cell affects the dynamics of the nucleus 

To assess nuclear dynamics independent of cell migration, we used micro-patterned fibronectin- coated substrates to 
confine cells to regions of defined geometry and size. Single cells were cultured on each patterned substrate and time 
lapse phase contrast imaging was carried out for about 8 hours (or till the cell underwent mitosis). Geometries with a 
variety of rotational symmetries circle, square, equilateral triangle and rectangles with aspect ratios 1:3 and 1:5 - but 
the same cell spreading area (1600/Ltm 2 ) -were fabricated and used to study effect of cell shape on the translational 
and rotational movement of the nucleus. Figure 1A-E shows color-coded intensity-profile images obtained by average- 
intensity projection of phase-contrast time lapse images for the above cases (Figure 1A-C are for rectangles with 
aspect ratio 1:1, 1:3 and 1:5 , Figure ID and E are for triangle and circle). The dark color or low intensity at the 
vertices of the triangular and rectangular patterns shows the formation of stable contacts in that region, a feature 
absent on the circular pads. Note that the cell adheres much more stably on the triangular pattern than on the circle 
or the square. 

The translational and rotational movements were measured from the time lapse images, and show convincingly the 
influence of cell geometry on nuclear dynamics. Figure IF displays typical trajectories of the nucleus on triangular 
and rectangular geometries of same area. On rectangles, presumably because of narrower confinement, the nucleus 
undergoes mainly translation while on triangles (as well as circles and squares) the nucleus both rotates and translates, 
as shown in Supplementary Movie 1-5. Since motion out of plane is negligible, we resolve the dynamics into two- 
dimensional translation and rotation in the XY plane. Translation is estimated through the instantaneous mean 
position (x,y) — { Xl 1 2 Vl , X2 + V2 ) of two nucleoli situated at roughly diametrically opposed points (xi,y\) and (x2,2/2)- 
The top inset to Figure IF shows a typical translation trajectory. Rotation is characterised by the coordinates of one 
nucleolus relative to this mean. A typical rotational track for the nucleus is shown in the bottom inset to Figure IF. 
Although rotation of the nucleus is not a normal feature of cell cycle, both translational and rotational movement of 
the nucleus during cell migration has been reported [11-18] for many cell types including NIH3T3 which was used in 
all our experiments. For completeness, we document such motion here as well. Figure SlA shows a representative 
DIC image of a monolayer of NIH3T3 cells cultured on glass bottom dishes. Time lapse images (Figure SIB) of three 
cells from this field of view are presented with arrows showing the position of the nucleolus. Rotation and translation 
tracks of the nucleus are plotted for these cells in Figure SIC and D respectively, showing large departures from its 
initial position and orientation. Finally, in order to demonstrate that migration is not the underlying cause of the 
rotation we observe, we confine cells by plating them onto fibronectin patterns of various well-defined geometries, 
allowing us to study the effect of cell geometry alone on nuclear dynamics. 

The fraction of rotating nuclei decreases significantly in geometries with large aspect ratio whereas on more sym- 
metric patterns namely equilateral triangles, squares and circles, the fraction is not significantly different (Figure 
2A). Except on rectangles, about 80% of cells showed at least 90° nuclear rotations in 8 hours. Nuclear circularity 
as a function of cell shape is altered with changes in aspect ratio (decreases from 0.9 to 0.7, Supplementary Fig- 
ure S2A) but not with changes in rotational symmetry of constraints. The instantaneous linear velocity decreases 
marginally from circle to square to triangle as well as with increase in aspect ratio of rectangle (1:1 - 0.23 //m/min 
and 1:5 0.20 ^m/min, see Figure 2B and Supplementary Figure S2B). The mean rotational velocity decreases from 
circle (2.2°/min) to square (2.0° /min) to triangle (1.6°/min) (Figure 2C and Supplementary Figure S2C, D and E) 
suggesting that rotation is sensitive to geometric constraints. However, the fraction of nuclei showing significant and 
systematic rotation is similar for these three shapes. To explore the possible role of myosin induced contractility in 
these phenomena, we turn now to the active hydrodynamic theory [19-23] of the cell interior. 

Active fluid with a central inclusion 

We show that cytoplasmic flows produced by acto- myosin contractility are the minimal explanation for the observed 
rotation of the nucleus. To this end, we turn to the theoretical framework of active hydrodynamics [26-29] . Contractile 
stresses carried by actomyosin, given an arrangement of filaments compatible with the cell shape imposed by the pads 
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and the presence of the nucleus as an internal obstacle, lead to organized flows that rotate the nucleus. More detailed 
propulsive elements, e.g., pushing by microtubules anchored onto the nuclear surface [7, 17], while possibly present 
in the cell, are not a necessary part of the mechanism. Since the cell in the experiment is stretched, its height is 
smaller than its dimensions in the plane. We can therefore model the cell as a quasi-two-dimensional film with the 
hydrodynamics being cut off at a scale proportional to the height. We also assume an axisymmetric cell, and ignore 
actomyosin treadmilling and the on-off kinetics of the motors. This highly simplified view of the cell still exhibits 
some key features of the dynamics found in the experiment. 

We now present the equations of active hydrodynamics [26-29]. The inner circular region represents the nucleus, 
which is taken to contain no active motor-filament complexes and is therefore modeled as a passive liquid drop of 
very high viscosity (rji) — in effect undeformable. The outer annular region is the cytoplasm, which contains active 
orientable filaments. The inner fluid-fluid interface, i.e., the boundary between cytoplasm and nucleus, has tangential 
stress continuity and tangential velocity continuity, and the outer surface, the contact line of cell with pad, has no 
slip. We assume the filaments preferentially lie parallel to any surface with which they are in contact. In particular, 
they therefore lie tangent to both the inner and the outer boundaries. 

The cytoplasmic medium is taken to consist of filaments suspended in the cytosol of viscosity rj <C r]i- We assume the 
filaments are in a state of well- formed local orientation whose manitude does not change so that it can be characterised 
completely by a unit vector or "director" field, n(r), [33] at postion r. Associated with the filaments is an active stress 
V4 7 n(r)n(r), where the parameter W is a measure of actomyosin activity, The concentration of filaments and myosin 
is assumed uniform. Fluid flow in the cytoplasm is described by the hydrodynamic velocity field v The equations of 
active hydrodynamics in steady state lead to a dynamic balance between shearing and relaxation of filaments [34], 
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with total stress tensor 
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Here F — J d d xf, f = (K/2)(Wn) 2 , is the elastic free energy for the director n, with the Frank elastic constant 
K, dfj — §ij — riiTij, Xkij — (l/2)(<5^7ij — S^jTii) + (X/2){8 r ^ i rij + S^-ni) is a flow-orientation coupling, and er° = 
(Vn)<9//<9(Vn) is the Ericksen stress. We will work with a completely symmetric stress, afj, built from cry, which 
will give the same velocity field, due to angular momentum conservation [35, 36]. 

The coefficient £ in (2) represents in a z-averaged sense the effects of confinement on the damping of velocities. It 
has two contributions: a viscous part ~ rj/h 2 arising because flows within an adhered cell of thickness h in the vertical 
z direction and no-slip at the base must in general have z-gradients on a scale h, and direct damping of flow through 
the kinetics of attachment and detachment of the cytoskeletal gel to the substrate. In our estimates below we retain 
only the viscous effect, so that C simply has the effect of screening the hydrodynamics at in-plane length-scales larger 
than h. Including attachment-detachment enhances £■ We use circular polar coordinates r, 4> in the plane. Since we 
assume axisymmetry, the radial velocity vanishes because incompressibility implies dv r /dr + v r /r — 0, and v r = at 
both the interfaces. For force balance in the region corresponding to the nucleus we have to solve the equation 
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where Q — rji/h 2 . The equation can be solved in terms of Bessel functions, with the constraint that has to be at 
r = 0. Continuity of tangential stress and velocity at the cytoplasm-nucleus interface gives the requisite number of 
boundary conditions. 



Force balance in the azimuthal direction reads 

d s 2n 
~dr a ^~ 
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Cv<p- (5) 



Expressing the diretor n = (n<^,n r ) = (cos 9, sin 9) the steady state equation (1) for the orientation field reads 

1 SF va, 

- -r- = — + A cos 29A r4j + ru = (A cos 29 - l)A r<j> (6) 
i 69 r 



5 



where A and $7 are the symmetric and antisymmetric parts of the velocity gradient tensor. 
Using (6), the r<p component of (3) can be recast as a first order differential equation for v^. 
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Thus, we have two first order equations, (7) and (5), and one second order equation (6) to solve, which we solve 
numerically 

A noticable and robust feature of the solution (inset to Fig. 5A) is the presence of a maximum in the magnitude 
of the velocity at some distance from the nucleus. This results from a combination of vanishing velocity at the outer 
boundary and the nuclear centre, and continuity of velocity and shear stress at the fluid-fluid interface. Note that 
our description does not include chiral effects, so that equivalent solutions with either sense of rotation are obtained. 

The competition between active stresses that promote flow and orientational relaxation, that inhibits it, is contained 
in the dimensionless combination a — W/(D ~ Wh 2 /K [27, 31]. Accurate estimates of parameters for our system 
are not easy to make. The cytoskeletal active stress, W, is generally argued to be in the range 50-1000 Pa [30]. Frank 
constants for actin nematics appear to be 2-20 pN [32], as in ordinary thermotropic nematics. The thickness h of 
the spread cell in our experiments is about 1/5 of the lateral extent. For a spread cell area of 1600/im 2 we therefore 
estimate h ~ 8/xm. Taken together, this leads to a ~ 200. However, if attachment-detachment contributions to £ 
are included, a will be lowered substantially. From the active hydrodynamic model we know [26] that the system is 
quiescent for small values of this parameter. However, a — 4.9, for which we present the results, is already sufficient to 
produce a spontaneous flow. Increasing a leads to increasingly complicated flows which we have only begun to explore 
[37]. We do not attempt a detailed comparison between the observed and the theoretical flow patterns. However, the 
conclusion about the maximum of the velocity being away from the nucleus rests purely upon the confining geometry, 
and we expect that the time and angle averaged velocity profile, measured from the experiment, will have a peak 
away from the nuclear boundary. 

For a > 1 i.e an unbounded, oriented active fluid, one expects [26] spontaneous velocity gradients of order 
In [27] and [31] it was shown that the presence of confinement on a scale h modifies the above conclusion giving 
a characteristic rate — F(j) [31] where F(x) — > 1 as a; — > oo and ~ x 2 forx — > 0, where i is the in-plane scale 
associated with observation. In our case, h is 1/5. Thus, the rotation rate should be of the order of Q.lW/rj. Using 
the arguments of [30] this estimate turns out to be of the order of a few degrees/min. This is reassuringly consistent 
with the magnitude obtained from the experiment. 

The two predictions we can make based on this simple model are that actomyosin is crucial for nuclear rotation, 
and that the angle and time averaged angular velocity will be maximum away from the nucleus. We perform a series 
of experiments to check these. In the next sections, we study the contribution of actomyosin contractility, a critical 
cytoplasmic regulator of nuclear prestress [38, 39], to the translational and rotational dynamics of the nucleus. 



We test the role of contractility on nuclear dynamics to validate the theoretical predictions based on active fluids 
with an inclusion. Actomyosin contractility was altered by treating cells fully spread on geometric patterns with low 
concentration of blebbistatin an inhibitor of the myosin II motor. To determine if the persistence in nuclear translation 
motion was dependent on contractility, the autocorrelation function (ACF), was plotted for control and blebbistatin 
treated cells (Figure 3A). Blebbistatin treated cells exhibited a decreased correlation time scale for translational 
motion (bottom inset to Figure 3A,) suggesting that actomyosin contractility is important for correlated translational 
movement of the nucleus. Next, the nuclear rotation angle as a function of time was calculated from the XY rotation 
trajectories. A typical plot of angle versus time for the nucleus on geometric pattern is shown in Supplementary 
Figure S3. On treatment with blebbistatin, the instantaneous angular velocity significantly decreases to 1.0 °/min 
when compared to control 1.6 °/min (top inset to Figure 3A and Supplementary Figure S4). To ascertain the effect of 
actomyosin contractility on the persistence of nuclear rotation we computed the auto-correlation of angular movement 
with time. Figure 3B shows plot of auto- correlation curve for control cells and cells treated with blebbistatin. Inset 
(below) to Figure 3B show that on perturbing actomyosin contractility, the persistence time of nuclear rotation 
decreases from 62 min in control to 33 min. In the next section, we study the role of actin flow patterns in regulating 
the nuclear dynamics. 



Role of actomyosin contractility on nuclear dynamics 
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Role of actin flow in driving nuclear dynamics 

Live cell fluorescence confocal imaging was carried out to simultaneously visualize actin flow dynamics and nuclear 
rotation on geometric patterns. Cells were transfected with lifeact-GFP to label actin in live condition (Figure 4). 
Time lapse confocal imaging of actin revealed a retrograde flow and its remodeling around the nucleus (Figure 4A) . 
To quantify the flow pattern (Supplementary Movie 6), we carried out particle image velocimetry (PIV) analysis using 
MatPIV. This revealed flow vectors tangential to the nuclear boundary with direction and magnitude correlated with 
that of the nuclear rotation as shown in Figure 4A and Supplementary Movie 7. Velocity field maps of actin flow were 
determined in small regions throughout the cell (Figure 4A, last panel and Supplementary Movie 8). A circulating 
flow, required to rotate the nucleus, is clearly seen (Figure 4A, middle panel and Supplementary Movie 6 and 7). 
Interestingly, upon blebbistatin treatment, inward flow of actin (Supplementary Movie 9-11), presumably driven by 
treadmilling, was not significantly affected. However, the azimuthal speed can be seen (Figure 4B) to decrease 
substantially, despite some scatter in the data. The circulation of flow around the nucleus was lost concurrent with 
the loss of nuclear rotation (Figure 4B, middle panel and Supplementary Movie 10). 

Further, we plot in Figure 5AandB, the angle averaged azimuthal velocity v^, with and without blebbistatin 
respectively, inferred from PIV as a function of radial distance from the centre of the nucleus. For comparison we 
also show the radial velocity, v r (Figure 5 C and D). Note that the graphs start from the edge of the nucleus. As 
predicted from the theory, the azimuthal velocity peaks away from the nuclear boundary in the control cells. In 
blebbistatin treated cells, by contrast, the velocity is 0, leading to the loss of nuclear rotation. However, v r is small 
in both cases, albeit with slightly larger fluctuations in the presence of blebbistatin. In addition, time lapse imaging 
of microtubules labeled with tau-EGFP show that the microtubule organizing centre (MTOC) undergoes translation 
dynamics while the nucleus exhibits both translational and rotational dynamics (Supplementary Figure S5). The 
orientation and arrangement of microtubules showed a cage like structure around the nucleus (Supplementary Movie 
12). This caging mechanism might help keep the nucleus relatively localized, thus enhancing the rotational effects of 
the torque generated by the actin flow. We also visualized the endoplasmic reticulum (ER) to assess its role in nuclear 
dynamics. Since ER is contiguous with the nuclear envelope, it could either stretch or undergo continuous remodeling 
as the nucleus rotates. Live cell imaging of ER, during nuclear rotation, showed dynamic remodeling suggesting a 
minor role for ER in nuclear rotation and reversals (Supplementary Figure S6). 

CONCLUSION 

Our results show that geometric constraints are critical in determining the rotational dynamics of the cell nucleus. 
While a number of components including cytoskeleton and motor proteins have been implicated to drive nuclear 
dynamics in migrating cells, our results on single cells confined to specific geometries suggest a role for actomyosin 
contractility. Square, circular, and triangular pads support mainly rotational motion, while long narrow geometries 
restrict it. The shape of the confining geometry further determines the magnitude of nuclear rotation; a relatively 
faster rotating nucleus is seen on circular pattern than on squares and triangles. We offer a simple theoretical 
explanation for the rotation in which the nucleus is modelled as a nearly rigid inclusion in the cytoplasm treated as a 
fluid containing filaments endowed with intrinsic stresses. The result is an angular velocity profile with a maximum 
at a radial position intermediate between the nucleus and the cell periphery, as observed in the experiments, and 
nonzero at the nuclear surface, corresponding to nuclear rotation. The predicted magnitude of the rotation rate 
based on plausible estimates of material parameters are also consistent with the measurements. That blebbistatin 
treatment greatly suppresses the flow lends support to our proposed mechanism. The question arises why nuclei are 
not universally observed to rotate in cells under normal conditions. At least two mechanisms could contribute to 
suppressing the generic instability that leads to circulating flows. One, in the absence of a rigid geometry, the cell 
boundary is free to change shape. This would disrupt the imposed boundary orientation of the filaments, and hence 
the orderly pattern of active stresses needed to drive a coherent flow. Two, the apical actin fibres [44], absent in 
square and circular geometries, present to some extent in triangular geometries, and very well formed in elongated 
geometries, bear down on and thus enhance the friction on the nucleus, suppressing its motion. 

Collectively, our results highlight the importance of both cell geometric constraints and actomyosin contractility 
in determining nuclear homeostatic balance. A number of experiments have shown that alterations in cell geometry 
affect gene expression programs [40] and cell cycle time [41], and lead to a switching of cell fates towards apoptosis 
or proliferation [42, 43]. We hope our work leads to a search for nuclear rotation in a wider range of systems and 
settings, whether such rotation has biologically significant consequences, and a deeper understanding of how the cell 
normally suppresses such effects. 
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FIGURES 




FIG. 1: Geometric constraints impinge on translational and rotational dynamics of the cell nucleus. (A-E) Colour -coded 
kymograph showing average intensity projection of time lapse phase contrast images for various shapes: rectangles of aspect 
ratio 1:1 (A), 1:3 (B) and 1:5 (C), triangle (D) and circle (E). (F) Typical XY trajectory of nucleus on triangle (red) and 
rectangle of aspect ratio 1:5. Inset: a typical translational (top) and rotational (below) trajectory of nucleus on triangular 
pattern. 
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FIG. 2: Quantification of the translational and rotational dynamics of the nucleus. (A) Fraction of nucleus which rotate on 
circles, squares, triangles, and rectangles of aspect ratio 1:3 and 1:5. (B) Instantaneous linear velocity of the nucleus on these 
patterns (C) Instantaneous angular velocity of the nucleus on circles, squares, and triangles. Error bars are SEM. 
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FIG. 3: Actin contractility affects persistence of both translational and rotational movement of the nucleus. (A) Mean 
autocorrelation curve for nuclear displacement in control (black) and blebbistatin (red) treated cells. Inset: Persistence time of 
translation for the two cases obtained by fitting the mean autocorrelation curves with single exponential decay function. (B) 
Mean autocorrelation curve for angular fluctuation for control (black) and blebbistatin (red) treated cells. Insets: Instantaneous 
angular velocity (top) and angular time scale (bottom) for the two cases. Error bars are SEM. 
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FIG. 4: Visualization of actin flow patterns during nuclear rotation in both control and blebbistatin treated cells plated on 
circular geometry. (A) Top panel: Tracks of two nucleoli (red and yellow) showing both translational and rotational dynamics. 
Scale bar = 20/im. Corresponding actin flow vectors (middle panel) and speed (last panel) was determined by particle image 
velocimetry (PIV) analysis using MatPIV for control (A) and blebbistatin (B) treated cells. The flow vectors have been scaled 
to 2 times for better visibility. Color code: 0.0 - 1.03 ^m/min. 
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FIG. 5: Azimuthal and Radial velocity of actin flow plated on circular geometry. Plot of and v r from velocity vectors of 
actin flow for control (A) and (C); and blebbistatin treated cells (B) and (D). Each color represents single time point for cells. 
Thick red curve is the mean of various such realizations (30 for control and 25 for blebbistatin treated cells) .Inset to 5A: A 
typical angular velocity vs. radial distance curve obtained by solving the equations (5), (6), and (7). 
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SUPPLEMENTARY INFORMATION 
Supplementary movies 

[Movies are available by e-mail to G.V. Shivashankar] 

Movie 1: Phase contrast images of the cell on circular geometry. Images were acquired every 5 minute. Scale bar= 
20 micron. 

Movie 2: Phase contrast images of the cell on square geometry. Images were acquired every 5 minute. Scale bar= 
20 micron. 

Movie 3: Phase contrast images of the cell on rectangular geometry of aspect ratio 1:3. Images were acquired every 
5 minute. Scale bar= 20 micron. 

Movie 4: Phase contrast images of the cell on rectangular geometry of aspect ratio 1:5. Images were acquired every 
5 minute. Scale bar= 20 micron. 

Movie 5: Phase contrast images of the cell on triangular geometry. Images were acquired every 5 minute. Scale 
bar= 20 micron. 

Movie 6: Actin flow pattern on circular geometry. 

Movie 7: Actin flow vectors on circular geometry. 

Movie 8: Actin velocity field maps on circular geometry. 

Movie 9: Actin flow pattern in cells treated with blebbistatin on circular geometry. 
Movie 10: Actin flow vectors in cells treated with blebbistatin on circular geometry. 
Movie 11: Actin velocity field maps in cells treated with blebbistatin on circular geometry. 
Movie 12: 3D image of microtubule constructed from high resolution confocal images. 



Supplementary figures 



15 




6 

? 3 
H. 

W 

X 
CD 

> 3 
-6 



0' 


6' 


12' 


18' 24' 30' 36' 42' 48' 


54' 


B 


B 


B 


B 






B 


B 




a 


a 


1 










MM 


B 


B 







































.*«* 


■#-*» 


>-* 




















* 


























1 




• 















































































































D 



E 


3 










m 





x 




CO 




> 


-3 




-6 





X axis (^iim) 



-6 -3 3 6 

X axis (fim) 



FIG. SI: Nuclear rotation in a monolayer of NIH3T3 cells. (A) Representative DIC images of NIH3T3 cells growing as a 
monolayer on glass bottom culture dishes. Scalebar = 100/im. (B) Time lapse images of single cells from different regions 
(cyan, red and green rectangles shown in (A)) of the monolayer. Time points are indicated at the top of each image. Plot of 
rotational (C) and translational (D) movement for the three cells shown in (B). 
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FIG. S2: (A) Nuclear circularity calculated from phase contrast image of the cell on various patterns. (B) Box plot of 
instantaneous linear velocity showing mean(small box), median (horizontal line within the box), 25th to 75th percentile (box 
range) and 5th to 95th percentile (vertical line) of nuclear movement on circles, squares, triangles, and rectangles of aspect 
ratio 1:3 and 1:5. (C) to (E) Histograms of instantaneous angular velocity (vertical box) and corresponding single exponential 
decay fit (red line) for circles, squares and triangles. The coresponding instantaneous angular velocities are quoted in the plot. 
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FIG. S3: (A) A typical angle versus time plot depicting various regimes (rotation, reversal and stalling) of rotational motion 
(red). Blue line joins the vertices obtained from Ramer DouglasPeucker algorithm for piecewise linear function which reduces 
the number of points from the curve depending tolerance. Linear fit to these regimes obtained from the above algorithm (green). 
Histogram of slopes (B) runtime (C) and runlength (D) for various regimes during rotation for multiple cells on triangular 
pattern (control-black) and blebbistatin treated cells (red). To further analyse nuclear rotation dynamics, the angle versus time 
was divided into different regimes: rotation, stalls and reversals. Linear fits were generated for the three respective regimes to 
compute the average angular velocity, runtime and runlength and their distributions were found to be similar for both control 
and blebbistatin treated cells. Although the typical nuclear rotational trajectory shows the three regimes, the ratio of total 
runtime in positive and negative directions is 1.8. 



18 




400 

Time (min) 



10 20 

56/5t (degree/min) 



50 

Time (min) 



BOO 



>> I 
o 




0)0.4 




w 








"§0.2 

to 


1 t 1)2 =1.0 degree/min 


"ra 

|o.o 


L 




100 




D 



400 

Time (min) 



800 



u —•—Control 
(U q j I — * — Blebbistatin 



| 0.1 

ra 

^ ,* 
2 o.o ittrlil* 




-200 200 

Angle (degree) 




50 

Time (min) 



100 



FIG. S4: Plot showing angle versus time curves obtained from nuclear rotation trajectory for control (A) and blebbistatin (B) 
treated cells on triangular geometry. (C) Histograms of instantaneous angular velocity (vertical box) and corresponding single 
exponential decay fit (red line) for blebbistatin treated cells. (D) Histogram of the rotation of the nucleus in both control and 
drug treated cases. Inset: Box plot of angle showing decreased width in the drug treated case. Angular autocorrelation curves 
for control (E) and Blebbistatin treated cells (F). The correlation timescale was determined by fitting individual curves to a 
single exponential decay. Each color represents different cells. 
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FIG. S5: Movement of microtubule organization centre (MTOC) labeled using tau-RFP. Tracks of MTOC (green), movement 
of two nucleoli (red and pink) overlaid on fluorescence and DIC images. 
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FIG. S6: Endoplasmic reticulum (ER) labeled using dsRED ER. Panel of images show reorganization of ER, near the nucleus, 
with time. 



